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Many strongly coupled fluids are known to share similar hydrodynamic transport properties. In 
this work we argue that this similarity could extend beyond hydrodynamics to transient dynamics 
through the presence of non-hydrodynamic modes. We review non-hydrodynamic modes in kinetic 
theory and gauge/gravity duality and discuss their signatures in trapped Fermi gases close to unitar- 
ity. Reanalyzing previously published experimental data, we find hints of non-hydrodynamic modes 
in cold Fermi gases in two and three dimensions. 


In the past decade, precision experiments of ultracold 
quantum gases at unitarity have enabled the study of 
transport phenomena in strongly coupled systems |lH6(. 
Usually, transport phenomena are encoded in hydrody¬ 
namic transport coefficients, such as the speed of sound, 
shear and bulk viscosities, heat conductivities, spin dif¬ 
fusion coefficients, etc. Presently, the speed of sound has 
been measured in a three dimensional unitary Fermi gas 
ii, and the shear viscosity has been constrained for 
both three dimensional and two dimensional Fermi gases 
0, 5] • The extremely low values of shear viscosity (when 
expressed in units of entropy) in particular suggest sim¬ 
ilarities between cold Fermi gases close to unitarity and 
very different systems such as hot quark gluon plasmas 
[7j, ;§j, high-temperature superconductors [§[ and strongly 
coupled fluids described by black holes via the AdS/CFT 
conjecture fiol ], all of which share similar transport be¬ 
havior. This apparent similarity in otherwise completely 
different physical systems suggests that these systems 
could be part of a broader class of so-called strongly in¬ 
teracting quantum fluids (SIQFs). 

It is conceivable that SIQFs share other properties be¬ 
sides their similar (hydrodynamic) transport behavior. 
This could be important because it could imply that it 
is possible to learn about one example of SIQFs (say 
high-temperature superconductors) through studying a 
different SIQF for which a particular trait is more easily 
accessible. In the present study we will investigate ultra¬ 
cold Fermi gases close to unitarity and argue that they 
exhibit properties similar to black hole SIQFs. 


One property that is quite remarkable about black hole 
SIQFs is that they do not seem to possess a description 
in terms of weakly coupled quasiparticles .11]- Instead, 
black holes can be characterized in terms of their ring- 
down s pec trum, similar to a glass struck (lightly) with 
a fork [l2| . Some of these quasinormal modes can be 
recognized to be the well-known hydrodynamic modes, 
i.e. sound and shear excitations. Others do not have 
an equivalent in (Navier-Stokes) hydrodynamics, and are 
thus non-hydrodynamic, but nevertheless affect trans¬ 
port properties (particularly on short time scales). 


If properties of SIQFs are universal, one would expect 
the presence of non-hydrodynamic modes in cold Fermi 


gases close to unitarity. This provides the motivation 
for searching for non-hydrodynamic modes in cold Fermi 
gases, both theoretically and experimentally. 

Transport in hydrodynamics and beyond: In 
order to understand the properties of non-hydrodynamic 
transport, let us briefly recall the properties of transport 
within hydrodynamics. To be specific, let us treat the 
case of a single component uncharged fluid in D spa¬ 
tial dimensions described by the Navier-Stokes equation 
(which clearly will not be applicable to Fermi gases below 
the superfluid phase transition temperature). The fluid is 
then characterized by the mass density p, the fluid veloc¬ 
ity u, the pressure P and the temperature T. Assuming 
hydrodynamic transport to be dominated by the shear 
viscosity coefficient 77 , we may set the bulk viscosity and 
heat conductivity to zero (this is a good assumption for 
Fermi gases close to unitarity, 0 ) . Finally allowing for 
a force term F (e.g. through a trapping potential) the 
fluid must obey the equations of motion 

d t p A di (put) = 0 , 

dt (frO j ) A dj (putUj A P&ij A rrtj ) pPi , 
d t £ A dj [uj (£ A P) A WijU) = pF -u, 

-Wij = nij , (1) 

with £ = -|- Rp an d (Ji - = (Q iU j A djUi — jjStjd ■ it). 

To close the system of equations, we adopt an ideal 
equation of state implying P = c 2 s (T)p , with c 2 (T) the 
temperature-dependent speed of sound squared (we as¬ 
sume an isothermal system). Additionally, we assume 
-p = const for simplicity. Considering small pertur¬ 
bations around some (possibly space-dependent) equi¬ 
librium configuration we have p = po( x ) + dp(t,x ), 
u = <5u(t,x), c 2 = Cq A 5c 2 (t). Hydrostatic equilibrium 
requires p^F = VPq. For perturbations around a con¬ 
stant density po, there are two solutions for the pertur¬ 
bations, namely the familiar hydrodynamic sound and 
shear modes. It is useful to characterize these modes in 
Fourier space, e.g. p ~ e —*wt+»k-x_ p or ^j ie sounc [ m ode 
(coupling perturbations dp, d ■ u) one finds the dispersion 
relation w = ±co|k| — i ^p Q CQ (l — -p), whereas for the 
shear mode (perturbations u transverse to k) one finds 

1 2 2 

w = —i ri p ) 0 . Thus, one expects to find density pertur- 
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bations which for fixed k behave as 


Wro(i.x) OC e ±-o|k|t+ l k.a : -ro t k 2 


( 2 ) 


where To = ^ (l — Eq. © was derived using hy¬ 
drodynamics, so its regime of validity is that of low k. 

The result for the sound mode perturbation in hydro¬ 
dynamics is to be contrasted with the result found for 
black hole SIQFs. Ref. 14] calculated the poles of the 
energy-momentum tensor correlator in the sound mode 
channel, finding the complete set of quasinormal modes 
for D = 3. For small k, the spectrum can be approxi¬ 
mated as 

nk 2 2 

U} h = ±c 0 |k| - ij— - , uinh,i ^ 27tT(±1.73 - 1.34*), 

4-to o 

( 3 ) 

and l = 27rTn(±l — i). The result for ujh is just 

that of a hydrodynamic dispersion relation for relativis¬ 
tic fluids in D = 3 (note that for Ref. 141 ^L. = _i_ 
and hence w n ?i,n>l = ^^(±1 — i) is also a consistent 
interpretation). The presence of the non-hydrodynamic 
modes u> n h, however, implies that density perturbations 
for fixed Ikl -C 1 behave as 


SpBH(t ,x) oc e ±-o|k| t +,k.x-iv + 0 nfc|n e ±flto "» t - r » t , 

n= 1 

( 4 ) 

with T = -Imw and amplitudes a n h , n >o for the non¬ 
hydrodynamic modes. Note that an infinite tower 
of quasinormal modes has also been found in non- 
relativistic systems for black hole SIQFs, cf. 0 ]. Also, 
non-hydrodynamic modes have been discussed before in 
the context of condensed matter systems in Ref. fid ]. 

Unless the amplitudes a n are extremely small or the 
damping rates T ra are very large, one could expect the 
non-hydrodynamic contributions to the evolution of den¬ 
sity perturbations to be experimentally observable. In 
the following, we will study the case of cold Fermi gases 
close to unitarity and investigate whether they exhibit 
non-hydrodynamic behavior similar to Eq. 0 - 

Non-hydrodynamic collective modes for a 
trapped Fermi gas: The analysis of collective modes 
changes if the equilibrium background density po( x ) is 
space-dependent, as is the case for experiments on Fermi 
gases in a trap. For an idealized harmonic trapping po¬ 
tential with base frequency w_l, the background density 
may then be written as 


logpo(x) = 


u]\{x 2 + y 2 + A 2 z 2 ) 
2cq 


( 5 ) 


where A = 0 for D = 2 and A « 1 for D = 3 (elongated 
trap). It is possible to classify all collective modes in 
this case by realizing that in the isothermal limit, Sc 2 (t ) 
is only a function of time, and hence the equations of 
motion place constraints on the space-dependence of V • 


hu. Thus, rather than analyzing perturbations in Fourier 
space, it is useful to expand perturbations in terms of a 
power series of x,y, and z. Restricting to inequivalent 
polynomials under rotations in the x,y plane, this leads 
to the ansatz 


— = coo(<)+cn(t)r cos(f>+C2o{t)r 2 +C22(t)r 2 cos 2 ^+... , 
Po 

( 6 ) 

with x = r cos </>, y = r sin q f>. Solving the equations of 
motion 0 then leads to the identification of three col¬ 
lective modes, corresponding to the time evolution of cn, 
Coo coupled with C 20 , and C 22 , respectively. These modes 
may be recognized to be the usual sloshing (“S”), breath¬ 
ing (“B”) and quadrupole mode (“Q”), respectively, and 
one finds the following dispersion relations in these chan¬ 
nels: 


w± 


UB ' h -±l 

/ 2 +l- 


(l- 1 ) 

W_L \ 

/ D 

Po 

V d) 


^A = ±V 2-^. (7) 

W_L Po 

All of these modes are driven versions of the sound mode 
excitations. One can explicitly verify that for an isother¬ 
mal gas with trapping potential 0 no shear mode per¬ 
turbations are excited. 

Unlike the non-hydrodynamic quasinormal modes en¬ 
countered in Eq. 0 , the hydrodynamic damping rates 
in Eq. 0 monotonically decrease in the ideal hydrody¬ 
namic limit 77 —> 0 . 

It is possible to try to extend the collective mode analy¬ 
sis for a trapped Fermi gas beyond Navier-Stokes hydro¬ 
dynamics. One option to do this is within the frame¬ 
work of second-order hydrodynamics 0 0. While 
second-order hydrodynamics is able to correctly capture 
higher order gradient corrections (e.g. 0(k 3 ) corrections 
in Eq. 0 ), its regime of validity is still uj -C 1. Thus 
one cannot expect second-order hydrodynamics to cor¬ 
rectly capture non-hydrodynamic modes. Nevertheless, 
one may try to analyze collective modes for a Fermi 
gas in a trap using second-order hydrodynamics to see 
if non-hydrodynamic modes emerge at least on a qual¬ 
itative level. This can be done easily since for linear 
response, second-order hydrodynamics basically replaces 
17 ij = ~ r l (T ij 77 ij + Tndt'Xij = —V a ij> with T n the re¬ 
laxation time for shear viscous stresses 0,0. Thus, 
one finds that in addition to the hydrodynamic modes 
in Eq. 0 . there is a non-hydrodynamic mode for the 
quadrupole mode oscillation in D = 2,3 and for the 
breathing mode oscillation in D = 3, given by 


^ B,nh ^ MQ,nh ^ 

^_L ’ CU± T n UJ ± ’ 


( 8 ) 


Another option, seemingly distinct from second-order 
hydrodynamics, is to consider kinetic theory in the form 
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of the Boltzmann equation. For kinetic theory with a 
simple collision term in the form of a relaxation time ap¬ 
proximation, the collective modes of a Fer mi g as in a har¬ 
monic trap have been analyzed before [3, 19|-21|. For a 
relaxation time tr in the collision term of the Boltzmann 
equation, one finds dispersion relations for the breathing 
and quadrupole modes given by 


iujq(wq — Aui\)tr + (2w^ — oJq) — 0 , D — 2,3 


iu B {wB - 4 u\)t r + ( ywi 




= 0, 0 = 3, (9) 


and = Au\ for 0 = 2. Eqns. d9j) have three solutions 
for ujq,ujb each, two of which agree with the hydrody¬ 
namic modes Eq. GD if ujj_tr = <C 1. However, in ad¬ 
dition to the hydrodynamic modes, Eqns. © contain one 
non-hydrodynamic mode for both the quadrupole mode 
and breathing mode in O = 3 (cf. [19j, |22|). In the strong 
coupling limit tr —> 0 , these non-hydrodynamic modes 
obey the dispersion relations © upon setting r n = tr. 
(Note that in second-order hydrodynamics 77 , 7 V are two 
different parameters, while in kinetic theory tr controls 
both the viscosity and the viscous stress relaxation time). 
The fact that the kinetic theory result matches that from 
second-order viscous hydrodynamics is expected since ki¬ 
netic theory is known to be a special case of second-order 
hydrodynamics in the strongly interacting tr —> 0 limit. 
However, one furthermore finds that in kinetic theory, 
non-hydrodynamic modes are present for all values of tr , 
describing a purely damped response in the quadrupole 
and D = 3 breathing mode channel. In particular, one 
finds from Eqns. © 


U B ,nh 

Ld± 


5 i 

StrLOjl 


(D = 3), 


^Q,nh 
LU± 


i 

2t r uj± ' 


( 10 ) 


in the limit of weak interactions tr — > 00 . Thus, while 
© is beyond the regime of validity in second-order hy¬ 
drodynamics, Eq. m is well within the regime of va¬ 
lidity for kinetic theory (weak interactions, well defined 
quasiparticles). 

Experimental evidence for non-hydrodynamic 
collective modes: Summarizing the theoretical status 
of non-hydrodynamic transport, it is known that non¬ 
hydrodynamic modes exist for some SIQFs such as black 
hole duals [lOj. Furthermore, for trapped Fermi gases it is 
known that non-hydrodynamic collective modes are fea¬ 
tured in descriptions beyond Navier-Stokes, but lie out¬ 
side the regime of validity of these descriptions for strong 
interactions [h| [23| . If the interactions are weak, and 
quasiparticle degrees of freedom are well defined, kinetic 
theory also predicts the presence of non-hydrodynamic 
modes, which are purely damped excitations m- 

In view of this, experimental input is crucial to hope to 
understand non-hydrodynamic transport for SIQFs. For 
this reason, we reanalyze the raw data for the collective 
oscillations observed in a D = 2 gas of 40 K atoms from 


Refs. @, [24| as well as raw data for the radial breath¬ 
ing mode oscillations observed in a D = 3 gas of 6 Li 
atoms close to unitarity from Refs. mm. To analyze 
the data for the time-evolution of the quadrupole (“Q”) 
and breathing (“B”) mode amplitudes, we use a fitting 
function of the form 


F(t) = a h cos(w h t)e 


—r ht 


■ a nh cos (w nh t)e 


-r„ 


( 11 ) 


where irrelevant phase shifts and offsets have been sup¬ 
pressed. Minimum y 2 fits to the data of Refs. !5j, I24pp26| 
are performed to extract frequencies w = Rew, damping 
rates F = -Imu as well as the amplitudes a for the hy¬ 
drodynamic (“h”) and non-hydrodynamic (“nh”) modes. 

We employed the Akaike information criterion (AIC) 
[27p to avoid overfitting the raw data with too many 
parameters. Using the AIC values extracted for both 
D = 2,3 data suggests that the currently available 
data is not sufficient to extract information about non¬ 
hydrodynamic frequencies. For this reason, we assume 
u n h = 0 in the following. 

For the case of D = 2, we analyze raw data by Vogt 
et al. in Ref. [B[ given at ^ = 0.45 with Tp = y/2Nco±, 
u >_l ~ 27t x 125 Hz for the case of N ~ 2000 40 AT atoms 
for four different values of the interaction strength pa¬ 
rameter ln(fcir'a). At each value of ln(fc^a), six different 
measurements of the time evolution are available and we 
obtain mean values and error estimates for the fit pa¬ 
rameters in Eq. EID from averaging over these six sets. 
For each of the ln(/cFa) values, we extract hydrodynamic 
mode frequencies and damping rates that are consistent 
with published values [51, |24j. 

In Fig. [TJ amplitudes of putative non-hydrodynamic 
modes extracted from the raw time-evolution data 
in Ref. [5] are shown, indicating a sizeable non¬ 
hydrodynamic mode component amplitude in the 
quadrupole excitations. We find that for a total of 240 
data-points, the AIC for the quadrupole mode decreases 
by 80 units when allowing for a non-hydrodynamic mode 
with a n h ^ 0 to be present in Eq. (fill) . The reduc¬ 
tion in the AIC value for the D = 2 breathing mode 
is similar, but the extracted non-hydrodynamic mode 
amplitude is consistent with zero (see Fig. ©. Our in¬ 
terpretation of this finding is that there is evidence for 
the presence of a non-hydrodynamic mode in the D = 2 
quadrupole mode data, while we find no evidence for such 
a mode in the D = 2 breathing mode data. The results 
for the extracted non-hydrodynamic mode damping rates 
are also shown in Fig. [l] Curiously, the extracted damp¬ 
ing rate ^ for In(kpa) ~ 1.84,2.8 is consistent with 
the kinetic theory analytic result m for an ideal gas 
in a harmonic trap in kinetic theory using the relation 
w±tr = 0.12 (l + ^ ln 2 (fcFo)) (cf. Refs.(l9l. 28|). 

For the case of D = 3, we analyze raw data by Ki- 
nast et al. in Ref. [ 25 } for the breathing mode oscil¬ 
lations of a gas of N ~ 2 x 10 5 6 Li atoms at a mag¬ 
netic field of B = 840 G (close to unitarity) for several 
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FIG. 1. Non-hydrodynamic amplitudes and damping rates for the Fermi gas in a trap in two dimensions (upper panels) 
and three dimensions (lower panels). Shown are results from our reanalysis of the raw time-evolution data from the D — 2 
experiment by Vogt et al. in Ref. [5] (upper panels), from the D = 3 experiment by Kinast et al. in Ref. [25| (lower panels), 
as well as the analytic results from kinetic theory with ideal gas equation ofstate in a harmonic trap © • 


different temperatures T/Tp where Tp = (37VA) 1 / 3 u;j_, 
A ~ 0.045, oj_ l ~ 2-7T x 1700 Hz. For each of the T/Tp 
values, we extract hydrodynamic mode frequencies and 
damping rates that are consistent with published values 
}25l |2fj. In Fig. © amplitudes and extracted damping 
rates of a putative non-hydrodynamic mode are shown. 
For a total of 1100 data points, the AIC decreases by 
80 units when allowing for a non-hydrodynamic mode in 
Eq. (fill) . Extracted amplitudes and damping rates seem 
for the most part uncorrelated as a function of T/Tp , 
but it is curious to note that for several values of T/Tp, 
the extracted damping rate is consistent with the kinetic 
theory result © using the relation lu±tr = 

Q. Our interpretation of these results is that the fits 
for the D = 3 breathing mode hint at, but do not pro¬ 
vide statistically significant evidence for, the presence of 
a non-hydrodynamic mode. 

Conclusions: We have argued that there could be 
a new class of strongly interacting systems (“SIQFs”) 
sharing similar transport properties beyond the now well- 
established value of shear viscosity over entropy density. 
If this is the case, then non-hydrodynamic quasinormal 


modes in SIQFs, theoretically well-established using the 
AdS/CFT framework, would imply the presence of non- 
hydrodynamic modes in all other SIQFs. We have stud¬ 
ied collective modes in trapped cold Fermi gases close to 
unitarity in D = 2,3 both theoretically and by reanalyz¬ 
ing experimental data. Our analysis hints at the presence 
of non-hydrodynamic modes in cold Fermi gases, but fur¬ 
ther experimental work would be needed to corroborate 
our findings. 

We believe that a study of non-hydrodynamic modes 
in ultracold quantum gases would open a new window 
into understanding transport properties of these systems 
and other SIQFs. This could help to found a new theory 
of “strongly interacting quantum matter” of importance 
in many subfields of physics. 
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